UNCLASSIFIED 


AD  NUMBER 

AD888796 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


EROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  T  &  E;  8  Nov  71.  Other 
requests  shall  be  referred  to  ONR,  Attn: 
Code  461,  arlington,  VA  22217. 


AUTHORITY 

ONR  Ltr .  ,  2  9  Aug  73 


THIS  PAGE  IS  UNCLASSIEIED 


Special  Report 


PROGRAM  BALLOC:  A  MIXED-WEAPON 
ALLOCATION  MODEL  AGAINST  MULTIPLE 
NON-OVERLAPPING  AREA  DEFENSES 


By:  R.  de  SOBRINO  B  J.  RIPPLE  N.  J.  LEMONS 

Distribution  limited  to  U.S.  Gov't,  agencies  only; 

•  2  NOV  1971*  Other  requests 

for  this  docuffl&at  aust  be  referred  to 
OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VIRGINIA  22217 

0^  ij/ikl 

CONTRACT  N00014-71-C-0015 


Reproduction  in  whole  or  in  part  is  permitted  tor  any  purpose  of  the  United  States  Government. 


STANFORD  RESEARCH  INSTITUTE 

Menlo  Park,  California  94025  •  U.S.A. 


Special  Report 


October  1971 


I 


5,- 

% 

a 


PROGRAM  BAUOC:  A  MIXED-WEAPON 
ALLOCATION  MODEL  AGAINST  MULTIPLE 
NON-OVERLAPPING  AREA  DEFENSES 


By:  R.  de  SOBRINO  B,  J.  RIPPLE  N.  J.  LEMONS 


Prepared  for: 

OFFICE  OF  NAVAL  RESEARCH 
DEPARTMENT  OF  THE  NAVY 
ARLINGTON,  VIRGINIA  22217 


CONTRACT  N00014-71-C-0015 


SRI  Project  8777 


Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United  States  Government, 


Approved  by: 

I,  J.  WEINSTEIN,  Director 
Systems  Evaluation  Department 

E.  J.  MOORE,  Executive  Director 
Engineering  Systems  Division 


Copy  No. 


32 


AUSTRACr 


Pri)^;iam  BALIjOC  is  a  muciel  designed  to  optimally  allocate  a 
mi xed-w.eapon  olfensive  lorce  over  a  set  ol  multiple-valued  targets 
deiended  by  multiple  nonintersecting  random  area  and  subtractive 
local  de tenses. 

The  program  will  ma.ximize  the  e.xpectod  damage  over  the  total 
target  set,  with  the  restriction  that  no  more  than  one  weapon 
type  be  allocated  to  any  randomly  deiended  island  or  target 
subso  t . 

This  report  presents  the  concept  on  which  the  model  is  based, 
and  describes  the  implementation  oi  this  concept  as  a  computer 
program.  It  also  contains  a  users  manual,  and  an  example  of  the 
use  ot  the  model,  employing  fictitious  data  to  avoid  the  neces¬ 
sity  of  security  classification. 


PREKACK 


This  rt'pori  is  submit  ti-d  undt-r  Contract  N00011-71-C-0015 
to  the  01 t ice  of  Naval  Research  under  the  task,  "Programming  a 
Strategic  Bomber  Allocation  Model."  This  report  contains  a 
description  of  the  model  as  well  as  a  users  guide  tor  running 
the  program. 

The  work  reported  here  has  been  performed  in  the  Systems 
Evaluation  Department.  Engineering  Systems  Division  of  Stanford 
Research  Institute,  under  the  supervision  of  Dr.  Ricardo  de 
Sobrino,  Staff  Scientist,  Mrs.  Barbara  J.  Ripple,  Mathematician, 
acted  as  project  leader,  and  Mrs.  Nancy  J.  Lemons,  Mathematician, 
did  the  programming. 
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I  INTRODUCTION  AND  PURPOSE 


During;  the  past  le\^  years,  models  used  to  evaluate  the 
e I tect iveness  ot  strategic  systems  have  been  highly  developed 
and  widely  used.  In  particular,  one  general  class  ol  such  models, 
which  will  be  called  here  the  lyP*-’"  become 

attractive  in  the  rapid  evaluation  of  relative  effectiveness  of 
alternative  systems.  In  addition,  these  models  are  useful  in 
investigating  the  sensitivity  of  the  results  to  the  values  of 
system  parameters,  and  to  the  interactions  of  offense  and  defense. 

SRI  is  currently  under  contract  with  OSD(SA)  to  develop 
techniques  to  improve  the  representation  of  airborne  strategic 
systems  in  effectiveness  models  of  the  aggregated  type.  As  a 
result  of  that  work  a  weapons-allocation  model  against  targets 
defended  by  random  type  area  defenses  with  certain  restrictions 
has  been  made  possible. 

The  purpose  of  this  report  is  to  describe  BALLOC,  a  weapons- 
allocation  model  programmed  for  the  Office  of  the  Chief  of  Naval 
Operations . 


For  a  more  complete  discussion,  see  Ref.  1.  References  are 
listed  at  the  end  of  the  report. 
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II  MODEL  DESCRIPTION 


A .  General 

I'he  allocation  model  described  here  has  some  unique  features 
that  are  not  available  In  currently  exlstliq^  models  of  the  ag^tre- 
gated  type,  such  as  CODE  50,  AEM,  and  similar  ones. 

The  model  assu.nes  a  number  of  corridors  (or  Islands)  de¬ 
fended  by  area  defenses  (fighters)  that  do  not  overlap.  It  also 
assumes  that  no  more  than  one  type  of  bomber  (more  specifically, 
payload  type)  will  be  used  within  each  island.  These  are  the 
main  limitations.  It  does  not  appear  easy  to  eliminate  the  first. 
Methods  to  remove  the  second  are  well  under  development,^  though 
this  can  only  be  done  at  the  expense  of  increased  running  time. 

Distinctive  features  of  this  model  are  as  follows: 

(1)  Bomber  penetration  probabilities  are  a  function 
of  the  attack  size  on  an  island  (or  corridor). 

(2)  They  also  vary  within  each  island  as  a  function 
of  the  penetration  distance;  in  fact,  they  may  be 
different  for  each  individual  target. 

(3)  They  are  accepted  as  inputs  in  tabular  or  functional 
form  and  may  be  obtained  from  detailed  simulations 
if  desired. 
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u. 


Problem  Formulation 


We  assume  here  that  a  set  ol  taryiets  is  defended  by  non¬ 
overlapping  multiple  random  fighter  area  detenses.  Since  the 
defense  coverages  are  nonoverlapping,  the  targets  are  in  fact 
divided  into  N  disjoint  subsets  (or  corridors).  Each  corridor 
IS  defended  by  a  known  random  fighter  defense.  In  addition, 
subtractive  local  de tense  may  also  exist  for  some  (or  all)  tar¬ 
gets  in  tl  corridor.  We  further  assume  that  each  corridor  will 
be  attacked  by  at  most  one  type  of  weapon. 

We  stiould  note  that  ttie  model  allocates  the  weapons,  not 
the  carriers.  This  self-imposed  restriction  is  needed  to  achieve 
acceptable  computation  times,  and  is  generally  assumed  in  models 
of  this  type.  Once  a  weapon  allocation  is  obtained,  weapons 
could,  in  the  real  case,  be  loaded  on  specific  carriers  that 
could  then  be  routed  to  specific  targets. 


th 


Let  B  be  the  number  of  type  j  weapons  allocated  to  the 

ij 

corridor.  The  penetration  probabilities  of  the  weapons  are 


given  inputs  in  tabular  form,  P  (B  ,  R) ,  where  R  is  the  distance 

iJ  iJ 

from  the  entry  point  to  each  island.  Thus,  and  in  general, 

penetration  probabi  ities  are  different  for  different  target.s. 

Then  for  a  given  value  of  B, .,  the  penetration  probability  of 

each  weapon  can  be  obtained  from  the  prescribed  penetration 

probability  function  P  (B  ),  which  in  turn  determines  the 

ij  ij 

damage  function  of  each  target;  the  optimal  allocation  of  B 

ij 

that  maximizes  the  total  expected  damage  function  of  the  island 
can  be  efficiently  generated  by  the  use  of  Lagrange  multipliers. 
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Let.  f  (B  )  be  the  total  expected  damage  lunction  to  the  i 
J-J  t.] 

island  caused  by  type  i  weapons.  We  note  that  1  is  a  lunction 

tj 

in  tabular  torm  since  it  corresponds  to  the  sum  oi  the  expected 

th 

duma^tes  ot  all  the  targets  in  the  i  c<>rridor  with  each  tar^^et 
having!  in  general  a  dilferent  damage  lunction,  because  the 
targets  either  belong  to  dillerent  exemplar  groups,  or  have 
different  penetration  probabilities  due  to  geographic  location. 

The  problem  may  be  expressed  mathematically  as  follows; 
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uhore  J  Is  the  total  number  ol  weapon  types,  and  U  equals  1 

th 

it  the  type  j  weapon  is  assigned  to  the  i  corridor,  and  equals 
0  o  t  he  rw  i  s  e , 

The  generalized  Lagrange  multiplier  method  is  readily 
applicable  to  this  problem.  Th.,  method  consists  of  two  parts: 

(1)  Fur  a  given  set  of  values  of  the  Lagrange  multipliers 

.  =  ( ',  ,•••,*,  ),  let  B'  =  (B.  )  be  the  weapon  alloca- 
1  J  i 

tion  matri.x  that  maximizes  the  Lagrangian: 


N  J 


N  J 


L(  .,B)  = 


E  E  -  -J  “lEE  Esj(v“ij) 

i=l  j  =  l  i=l  j  =  l 


For  the  i  corridor,  let  j.  be  the  weapon  type 


that  maximizes 


L.  B.®\  , 

\  J  tj  /  ,  i.e. , 

L  /  .  ,B  ■’\  =  max  L  /  ,  ,B 

J  ijj  ij\  j  ij/ 


LetU'’  =  lifJ=j 

ij  i 


=  0  otherwise. 


i.e.,  assign  weapons  of  type  to  the  i 
corridor  and  let 


N 

B*  =  U  '  B  ' 

j  Z _ I  1.1 

x=l 


Tlien  U  '  and  B  '  are  optimal  for  the  assit^nment 

i.l  1.1 

pi'oblem,  \M  t h  B*  in  place  of  B  for  the  total 
number  ot  type  j  \^eapons. 


(2)  Using  an  existing  iterative  procedure,  the  values 
of  .  may  be  systematically  adjusted  so  that  the 
\xeapon  resource  vector  (B*)  converges  to  its  desired 


values  (B  )  . 

J 


III  MODEL  IMPLEMENTATION 


A .  I ntruduction 

The  proj^ram  consists  basically  ol  three  steps:  (1)  lindinn 
the  damage  functions  for  each  attacker  type  on  each  corridor, 

(2)  optimally  allocating  the  offense  force  on  the  corridors,  and 

(3)  adjusting  the  optimal  allocation  to  be  identical  with  the 
desired  allocation  when  the  latter  is  not  exactly  obtained.  This 
last  step  will  be  denoted  by  "closing." 

B.  Generating  the  Corridor  Damage  Curves 
1 ,  General 

Because  the  probability  of  bomber  penetration  to  each 
target  (and  the  corresponding  damage  function)  depends  on  the 
total  attack  size  to  a  corridor,  duo  to  the  assumption  of  random 
area  defenses,  the  following  steps  are  necessary  to  obtain  each 
point  on  the  damage  curve  of  a  given  weapon  type  on  a  given 
corridor: 

(1)  For  each  desired  attack  size,  find  the  penetra¬ 
tion  probability  corresponding  to  that  attack 
si  ze , 

(2)  Find  the  individual  damage  functions  for  each 
target  given  the  penetration  probability. 
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(3)  Usln({  LuKrauK*^'  multipliers,  tind  the  optimal 
alloeation  tor  the  desirt*d  attack  size. 

Kach  ot  thi'Se  thni  stips  is  discussed  in  more  detail 

be lou . 

The  number  ot  points  on  each  curve,  as  mil  as  the 
maximum  attack  si/*-  tor  lach  curve,  are  user  Inputs.  Due  to 
computer  storage  requ  i  remi  nt  s ,  hovu  ver,  a  maximum  ot  20  points 
per  curve  cannot  be  presently  exceeded. 

2 .  Finding  Penetration  Probabilities 

To  find  the  probability  of  a  weapon  penetrating  to  a 
target,  t,  the  great-circle  distance,  d,  from  the  corridor  entrance 
to  that  target  is  found  using: 

c  =  cos  [sin(LA  )  sin(l.^  )  +  cos(LA  )  cos(LA  )  cos(L0  -  LO  )] 
pi  T  pi  T  T  pi 

e  =  cos  ^  [sin(LA  )  sin(LA  )  -f-  cos(LA  )  cos(LA  )  cos{LD  -  LO 
pi  p2  pi  p2  p2  p 

B  =  cos  ^  IslnCLA  )  -  cos(c)  sin(LA  )]  'lsin(c)  cos(LA  )]  - 
T  pi  pi 

cos  ^  [sin(LA  )  -  cos(e)  sin(LA  .)]  'lsin(e)  cos(LA  )] 
p2  pi  pi 

d  =  IMOO  sin^[sin(c)sin(B)]  (1) 

where 

LA^,  LO^  =  Latitude  and  Longitude  of  target  t 

LA  ,  LO  ,  LA  ,  LO  -  Latitudes  and  Longitudes  of  cor- 
pT  pi’  p2  p2 

ridor  entry  boundaries 

d  =  Great-circle  distance  between  tar¬ 
get  t  and  corridor  mtry  line  de¬ 
fined  by  pi  and  p2. 
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^  Given  this  distance  and  the  attack  size  ot  interest, 

a  table  lookup  procedure  is  used  tt>  obtain  tlu'  penetration 
probability.  These  tables  are  inputs  that  must  be  generated 
e.xternal  ly .  "  Linea,'  i  nterpolat  ii>n  is  used  ti.)  obtain  the  desired 
attack  size  trom  the  tables. 


3 .  Target  Uamat;e  Functions 

The  square- root-law  damage  lunction  is  used  lor  each 
target.  For  a  locally  undetended  target  with  penetration  proba¬ 
bility  equal  to  one,  this  is  ol  the  following  form: 


where 

V  =  Target  value 
i 

a.  ,  =  Number  of  type  J  weapons  arriving  at  target  i 
J  >  ^ 

k  =  A  damage  constant  that  depends  on  the  target's 

hardness  and  size  and  the  weapon's  yield  and  CEP 


b  =  Expected  damage, 
i 


Given  a  penetration  probability  (P  )  for  attacker  i 

ij 


on  target  i,  a  new  value  ot  k,  k^  is  found  such  that  | 


I'hiM  k'  value  then  has  the  effect  of  "stretching"  the 
undefendeil  tlamaKc  function  to  aIlo«  lt>r  the  effects  of  the 
penetratiott  probability  throu^jh  the  area  defenses,  and  the  ex¬ 
pected  dama>;e  is  found  by  usinn  this  k'  value  in  Kq.  (2). 

It  tfierc  IS  a  local  defense  at  target  i,  the 

term  in  Kq.  (2)  is  replaced  by  yj .  .  "  Y,  ,  Mhere 

J  »  i  .1 » 1 

Y  =  I  R  k  P  (4) 

.),i  ,I,i  J  J,i 

and 

Y  =  Equivalent  price  of  the  target  i  for  weapon  type  j 
J  .  i 

I  =  Interceptors  at  target  i  effective  against  weapon  j 
J  ,  1 

R  =  Interceptor  kill  probability  against  weapon  type  j 
.i 

k  =  An  input  constant  associated  with  the  leakage 
through  the  area  defenses"' 

P  =  Penetration  probability  through  the  area  defenses, 

J  >  ^ 

a  function  of  the  target  location  and  attack  size. 

4 .  Find  the  Desired  Attack  Size 

Given  the  individual  target  damage  functions  at  a 
corridor,  it  is  then  possible  to  obtain  the  attack  size  of 
interest.  This  is  done  using  Lagrange  multipliers  in  an  itera¬ 
tive  technique.  Starting  with  a  very  large  value  of  a  single 
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scalar  multiplier  aiul  a  very  small  one,  attack  levels  that 

bracket  the  desired  attack  size  are  obtained.  A  neM  value  lor 
the  multiplier  is  found  by  takiUK  the  avera^je  of  the  tMo  previous 
ones.  The  neM  allocation  is  obtained  and  replaces  the  old  point 
on  the  same  side  of  the  desired  allocation.  This  avera({in({ 
procedure  continues  until  either  the  desired  point  is  found  or 
the  values  of  the  multipliers  do  not  coalesce  into  a  sin^jle 
value  resulting  in  the  desired  number  of  weapons,  in  which  case 
the  solution  is  in  a  "nap,"  When  this  occurs,  the  point  closest 
to  the  desired  point  is  used. 

C .  Optimal  Weapon  Allocation 

To  optimally  allocate  the  offensive  force,  (l.e.,  maximize 
the  expected  damage  destroyed  over  all  the  corridors)  the  program 
again  utilizes  the  Lagrange  multiplier  method.^  In  this  case, 
the  Lagrange  multiplier  is  a  vector  in  j  space  where  j  is  the 
number  of  weapon  types. 

Given  the  multiplier  vector,  the  problem  of  finding  the 
optimum  force  allocation  is  quite  simple,  since  the  damage 
functions  are  in  tabular  form  and  the  corridors  do  not  overlap. 
Each  corridor  can  be  treated  independently,  and  the  "best"  weapon 
type  (and  the  number  of  weapons  of  this  type)  to  use  on  the 
corridor  is  found  by  searching  each  damage  table  to  find  the 
point  that  maximizes  F(A)  -  A  =  L  for  payload  type  j,  and 
then  choosing  the  type  with  the  largest  L,. 
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Aftcr  each  allocation,  the  multiplier  vector  is  adjusted 
until  either  the  desired  attack  size  (to  Mlthin  the  granularity 
of  the  damage  table)  is  obtained,  or,  the  change  in  the  values 
ot  the  multipliers  is  very  small,  or  a  maximum  number  of  Itera- 
tlon  is  reached.  If  the  desired  attack  size  is  obtained,  the 
run  is  complete;  otherwise,  a  closing  procedure  is  required  to 
adjust  the  optimal  attacks  to  mutch  the  desired  one. 


To  adjust  the  multiplier  vector,  the  folloi^ing  standard 

algorithm  is  usod.»  Along  with  the  X  associated  with  each 

J 

attacking  payload  type  is  a  AX  .  It  is  the  that  is  adjusted 

J  J 

to  find  the  new  \  .  If  ^X '  is  going  in  the  right  direction  but 
J  J 

the  obtained  allocation  is  still  not  at  the  desired  value,  then 

Ax '  =  c  ^X',  where  c  is  an  input  constant  >  1.0.  If  in  the  last 

J  1  j  I 

iteration  the  desired  allocation  was  crossed,  the  AX^  =  c  ^X  , 

,i  2  j ' 

where  c  is  also  an  input  constant,  negative,  and  <  l.O.  Also, 

2 

a  Ax  term  is  applied  to  all  X  's.  This  term  is  modified  at 

each  iteration  by  the  same  constants  as  the  Ax^'s,  depending  on 

whether  Z  X^  W  is  greater  or  less  than  Z  X.W  (W  =  desired 
J  J  D  j  j  A  D 

allocation;  W  s  obtained  allocation).  Thus, 

A 


X'  =  (I  +  AX')  (1  +  Ax)  X 
J  i  j 


The  values  of 


c 

1 


and  c  have  been  obtained  heuristically . 


*  This  algorithm  was  originated  by  Dr.  H.  Everett  of  Lambda  Corp. 
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u. 


riosiiin  Pro ce dure 


The  cluslii^;  procedure  is  a  heuristic  method  used  to  slightly 
change  an  obtained  optimal  allocation  to  reach  the  desired  one 
'^hi'it  thi  latter  lies  in  a  "gap"  and  is  tiot  obtainable  by  the 
Lagrange  multiplier  method. 

In  Program  BALLOC,  procedures  are  needed;  the  lirst 

aitjusts  the  allocation  among  the  corridors  until  the  desired 
allocation  is  obtained.  The  second  procedure  tlien  adjusts  the 
individual  corridor  allocations  among  the  targets  within  that 
Ci^rr  idor . 


The  first  step  in  performing  the  first  closing  is  to  select 
the  "best"  allocation  to  close  on.®  This  is  done  by  finding  the 
allocation  that  minimizes: 


(5) 


\^here 


W 

A 

j 


=  Obtained  number  of  type  j  weapons 


W  =  Desired  number  of  type  j  weapons. 

J 


The  routine  imposes  the  restriction  that  each  corridor 

must  use  the  weapon  type  assigned  in  the  optimal  allocation; 

thus,  for  each  type-J  weapon,  if  Wa  ^  "V,  ,  the  routine  selects 

J  j 

from  the  corridors  that  are  assigned  this  type  the  one  receiving 
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tho  loast  damamo  per  attacker,  and  reduces  the  allocation  to 

this  corridor.  Similarly,  if  W,  <  ,  the  corridi>r  that  Mill 

j 

receive  the  most  damage  per  attacker  is  selected  to  receive  the 
additional  Meapons. 

If  a  Meapon  type  Mas  not  allocated  at  all,  tlu'  routine 
selects  from  the  unat tacked  corridors  the  one  that  Mould  receive 
the  most  damage  and  allocates  the  Meapons  to  that  coriidor. 

After  the  exact  allocation  tti  each  corridor  is  obtained,  tlie 
procedures  described  previously  in  Section  B-d  are  used  to 
allocate  the  Meapons  among  the  targets  Mi  thin  the  corridor. 

If  the  desired  allocation  to  that  corridor  is  not  exactly  ob¬ 
tained,  an  algorithm  similar  to  that  used  in  the  first  closing 
is  used  for  slight  adjustments. 
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IV  PROGRAM  UESCRIPTION 


A. 


I ai roduc  t ion 


BALLOC  iri  written  i  tt  Kortran  IV  lor  use  on  the  SRI  CDV 
6-100  computer  as  well  as  the  CNA  COC  BGOO.  Th.e  cnri-  storaei- 
reciuirement  is  approximately  50,000  words.  The  running  time 
can  be  broken  into  two  parts--the  time  required  tt>  generate 
ttie  damage  curves  tor  each  corridor,  and,  given  the  curves,  the 
time  required  to  obtain  a  given  allocation.  To  generate  a  curve 
with  20  points  tor  one  weapon  type  on  one  corridor  requires 
approximately  10  seconds  on  the  CDC  6-100.  The  running  time  to 
obtain  the  desired  weapon  allocation  depends  on  the  number  ol 
corridors  and  the  number  of  weapon  types.  For  two  weapon  types 
and  nine  corridors,  the  time  is  approximately  5  seconds. 


Section  IV-B  contains  flow  charts  of  the  main  program  and 
of  the  more  complicated  subroutines.  Section  IV-C  gives  a 
description  ot  the  user  input  to  the  program,  along  with  typical 
sample  values.  In  Section  IV-D  are  several  figures  showing 
various  representative  sections  of  the  output  of  this  program. 

The  target  values  and  damage  functions  used  to  calculate 
the  e.xamples  shown  are  tictitious  to  avoid  security  classification. 


This  is  an  average  value  assuming  40  targets  per  corridor  and 
obtained  from  a  total  time  of  about  180  sec  for  400  targets  lo¬ 
cated  in  9  corridors. 


B. 


Pro^;ratn  F I ou 


Fi^jure  1  is  a  flo\k  chart  of  the  main  program  flow.  Figure 
2  shows  the  allocation  subroutine,  and  Figure  3  the  closing 
procedure . 


C .  Data  Input 


1 . 


Card  ID 
1 
2 

3 


4 

5 

6 


8 

9 

10 


Data-Card  Deck  Set-up 

The  following  cards  are  required  to  run  the  program; 


Number  of  Cards 


Descri ption 


1  Variable  size  card 

1  Number  of  points  on  damage  curves 

for  each  weapon  type 

1  Maximum  number  of  weapons  of  each 

type  to  be  used  in  computing 
damage  curves 

1  \  adjustment  card 

A  set  for  each  "k"  values  for  each  exemplar  group 

weapon  type 


One  for  each  local-  Local-defense  card 
ly  defended  target 

1  Numbei’  of  targets  in  each  corridor 

One  for  each  Target-description  card* 

target 

One  for  each  Corridor  entry  points 

corridor 


1 


Interceptor  constants 


*  These  cards  are  ordered  so  that  all  Corridor  -1  cards  are  first, 
followed  by  all  Corridor-2  cards,  etc. 
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,  BALLOC 


BEAD 

INPUTS 


FIND  DlSTANCfc  FROM 
entry  points  to  targets 
Find  equivalent 
target  prices 


HE  AO 

penetration 

probability 

tables 


Find  damage  function 
for  each  weapon 
TYPE  ON  each 
CORRIDCR 


optional 

PRINTOUT 


READ  IN  A 
WEAPON 
ALLOCATION 


initialize 

LAGRANGE 

multiplier 

value 


■7=^3—. 

1  allocation 

update 

LAST 

MULTIPLIER 

ITERATION  ^ 

values 

OPTIONAL 

PRINTOUT 


PERFORM 

CLOSING 


PRINT 

RESULTS 


FIGURE  1  BALLOC  FLOW  CHART 


TA-8777-1 


ALOCAT 


INITIALIZATION 


S6T  CORRIDOR 
lO  TO  t 


FLAG  This 
WEAPON  TYPE 
FOR  this 
CORRIDOR 


SET  WEAPON 
TYPE  ID  TO  1 


FIND  POINT  ON  CURVE  FOR 
THIS  WEAPON  TYPE  ON  THIS 
CORRIDOR  with  largest 
LAGRANGIAN 


/  IS  N. 

This  n. 

LAGRANGIAN  > 
THE  ONE  FOR  ANT 
PREVIOUS  WEAPOI^ 
\  TYPE 


INCREASE 
WEAPON  TYPE 
COUNTER 


MORE  ' 
WEAPON 
TYPES 


'  MORE  ^ 
CORRIDORS 


INCREASE 

CORRIDOR 

COUNTER 


OPTIONAL 

PRINTOUT 


RETURN 


TA-8777-2 


FIGURE  2  ALLOCATION  SUBROUTINE 


Card  ID  Number  of  Cards 


Description _ 

1 1  I  Interceptor  reliability  attalnst  each 

weapon  type 

12*  1  Attack-size  table 

13*  As  many  as  needed  Distance  from  entry-line  table 

14*  As  many  as  needed  Penetration  probabilities  as  a 

function  of  the  attack  size  and 
distance 

15  One  for  each  case  Desired-allocation  card 

to  be  run 

2 .  Input-Card  Formats 

Figures  4  through  18  describe  the  input  card  formats. 

The  first  row  contains  the  number  type  (I  =  right  iustified  integer, 
R  =  a  real  number,  and  A  is  an  alpha  format).  The  second  row 
contains  the  field  on  which  the  number  is  to  appear.  The  third 
row  describes  the  Input,  and  the  fourth  row  gives  sample  values. 

Figure  4  contains  the  size  limitations  for  the  various 
parameters  as  well  as  giving  three  print  options.  If  the  Option 
I  field  is  not  blank  or  0,  damage  curves  for  each  weapon  type 
on  each  island  will  be  plotted.  Option  2  will  print  out  the 
allocation  to  the  corridors  for  each  iteration,  unless  the  field 
is  blank  or  0.  Similarly,  a  blank  or  0  for  Option  3  will  sup¬ 
press  the  printing  of  the  allocation  summary  for  each  cycle. 

Figure  7  contains  the  constants  for  increasing  or  de¬ 
creasing  the  Lagrange  multiplier  vector  to  converge  on  the  desired 
solution,  as  well  as  the  cutoff  point  to  stop  the  iterations, 

* 

These  are  inserted  In  the  above  order  for  each  corridor. 


Each  target  is  assigtu-d  to  an  exemplar  group,  hoaever, 
there  may  be  up  to  -100  groups  so  that  each  target  can  be  treated 
Independently.  The  damage  constant  ("k")  (Flgun-  8)  is  given  for 
each  lAeapon  type  and  each  exemplar  group. 

I'he  target  cards  (figure  ll)  are  ordered  st>  that  all  of 
the  targets  in  Corridor  I  are  first,  followed  by  Corridor  2,  etc. 
Figure  10  gives  the  number  of  targets  in  each  corridor.  The  last 
corridor  contains  all  of  the  targets  not  covered  by  any  area 
defense . 


For  each  city  Vkith  local  defense,  a  card  (Figure  9) 
shoviing  the  number  of  interceptors  against  each  vieapon  type  is 
used.  Also,  a  constant  determined  by  the  leakage  through  the 
area  defense  Is  needed  for  each  Vkeapon  type  (Figure  13).  Inter¬ 
ceptor  kill  probabilities  are  given  against  each  weapon  type 
(Figure  14),  and  to  obtain  the  equivalent  intercepts  of  weapon  J 
against  target  i  the  program  uses 


E 

IJ 


1  R  K  - 

ij  j  j  P 


P 


where 


I  =  Interceptors  at  target  I  against  weapon  type  j 
ij 

R  =  Interceptor  kill  probability  against  weapon  type  j 


K  =  An  input  constant  associated  with  the  leakage  through 

j 

3 


the  area  defenses. 


P  ,  t  hi-  piiutjatu»n  piobabi  1  1 1 V  Id  a  tarnct,  is 

P 

a  uiiu  tiDM  "t  tDtal  attack  si/c  tn  tbi-  cui  i  idor  and  the 

di>taiu'v  Dt  the  target  t  tdiii  the  corridor  entrance.  Figures  16, 

17,  anil  Ik  lontain  the  tabl«.‘s  dt-scribinn  these  t  iiin.' t  i  ons . 

D .  Sample  Output 

Figures  19  throui;h  23  are  samples  of  the  output  available  from 
BALLOC.  Fixture  19  contains  the  pro|»iam  inputs.  Figure  20  sHomk 
the  plots  ot  the  individual  damage  functions  tt>r  each  weapon  type 
on  eacii  corridt>r.  This  printout  i-  optional,  as  is  Figure  21, 
which  contains  the  output  available  after  each  iteration.  Figures 


FIGURE  5  DAMAGE-CURVE-FOINT  LIMITS 


up  to  5  weapon  types 


up  to  400  exemplar  group* 


FIGURE  13  INTERCEPTOR  CONSTANTS 


figure  15  ATTACK-SIZE  TABLE 


up  to  12  di  s  laiicus 


FIGURE  17  PENETRATION-PROBABILITY  TABLE 
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^IKal  allocation  TC  TAROetS  IN  CCRBIOCR  2 
T^IS  CORRIDOR  HAS  ATTACNCO  BY  WC«ROK  1 


TCTal  ATTACHfRS  TO  CORRICCR  ■  21 

TCTal  0AH46C  TO  CORRIDOR  •  in9.9SSl 
TCTAL  VALUE  OR  rORRiOOR  •  llSOO.OO 


T«R«CT 

ATTACK 

camaof 

3 

n 

0.0000 

2 

26,A2T? 

?1 

A 

120.2T52 

T'' 

T 

109.1215 

1? 

1 

13.5355 

SA 

2 

26.0A65 

S5 

0.0000 

A2 

r 

0.0000 

AA 

n 

0.0000 

T1 

1 

12.5532 

AT 

n 

0.0000 

inn 

n 

0.0000 

l«A 

n 

0.0000 

120 

f> 

0.0000 

123 

rt 

o.ooon 

129 

0.0000 

1  32 

0 

0.0000 

1  A9 

r 

0.0000 

ITS 

0.0000 

2n0 

n 

0.0000 

2''T 

A 

0.0000 

P'^A 

A 

0.0000 

212 

n 

o.oooo 

236 

n 

0.0000 

2A5 

n 

0.0000 

2a9 

ft 

0.0000 

250 

n 

0.0000 

2T9 

n 

0.0000 

29A 

n 

0.0000 

333 

r 

0.0000 

3An 

0 

o.oooo 

359 

0 

0.0000 

371 

0 

0.0000 

3A3 

0 

0.0000 

3P9 

n 

0.0000 

39A 

0 

0.0000 

399 

r 

0.0000 

FIGURE  22  FINAL  CORRIDOR  ALLOCATION 
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so»*»'*Bv  or  riN-»i  riLOCAtiON 


CCRBICOP 

HEAP  TTP 

attack 

Damage 

VALUE 

1 

n 

0 

0*00 

3*00 

? 

1 

?\ 

3n9,96 

11300 

3 

n 

0 

0  *00 

16920 

4 

f\ 

0 

0*00 

23760 

5 

f) 

0 

0*00 

10360 

t 

r 

0 

0*00 

7260 

7 

r 

0 

0*00 

6200 

e 

r 

0 

0*00 

6360 

9 

? 

37 

4497*06 

23200 

hEAP  TVP 

attack 

DAMAGE 

1 

3«9.9591 

? 

?7 

4497, 0R2G 

TCTal  A7TACKE«S  ■  5« 

TCTAL  damage  •  4en7,,74l9 

TCTAL  value  ■  110200*0000 


FIGURE  23  SUMMARY 


P«0GO4i-  PALIOC  nKP»JT.OUTPuT*TAP€^«TMPUT*TAPr<k«OUTPUT» 

C 

C0MMCK/rrnAp/LAPrA*t (SOOtS) #v(S00) trPTTLMi^rtOtS) * ifirst*tlasTi 
UTTAfX.OAMAfiFtATT (500J  *OAM<500) tNlNTfP (SOft J 
COMMrK/r»Ll  nc/w0(%)  tLAMPOA  (5)  *  ITYPF  ( 1H)  ^awf  AP(15)  tCnAMd^i  • 

UTTCBV  ( tTAMCPV  ( lS*S*2n  •LtMMAlt  (15*^)  *!NTE»V(1S)  # 

2T0TAI  A  (^)  *TrTAL0(5)  t  TOT  ALL  (S)  •VALlIP  (1^) 

COMMC^yC9€«;T/Tv<J•TOTA0^5)  *!TyPE8US>  .PwFAPPd^)  .EnAMedS)  • 

1  INTBVP  » 10)  *1  Ami>aB(S)  tlTEPPtlTEP 
C 

OlMfKSTON  NTAOr,  ( 1^)  .TO  1ST  (SOO)  #  I^PINT  r  15)  *  AP  ( 5)  •OIST  U  P )  * 
IPP(1?*^«5J#A(2!i5).NP(?1«5J  tKAY (5oo,5) t IPPfWT (3) *Wf APmy (15) • 
2inTST<B00J.X{2l)*Y{2n  t0LAP0A(5) * SAM {500*5) *LAT (500) • 

3L0N6(5PO)  .CLAT  (15*?)  •CLOMOdSta)  tKAYl  (400*5)  *PEL(5)  * 

410(50'’)  tNTKr»(«i)  *NINC»P(5)  tINCP(5)  .NWFPMXf5)  *0(5)  *1NTEP(5) 

c 

PPAL  t'AVtKvLAMOAl  t  L  4M0  A?  ^L  AMOA  tL  Amrd  A  •  L  AMma  V  ,  L  TQT  AL  *L  AmQAB  « 

lLATtLONO*KAvl 

C 

PE  AO  10'’*  NrOBin,NATTYPfNAPtNDIST,lTFPMX* ( TPPTKT ( ! ) * I«1 *3) tNEXEMP* 
INLOEF 

IF  (N ATTYP.OT.KroPtO)  GO  TO  650 
PEAO  )00*  (MNCP  (IW)  *  !«■)  tNATTYP) 

PEAD  loo*  (KMEPMX (Itt) •lM«i«NATTYP) 

00  ?nl  I*(«1,NATTvP 
nTmcbp  (TM)*MNrR  (1 4)* I 

20  3  iMCP  ( Tw)  ■  (Float  (N«iFPMx  ( iw) )  /nimcp (  tw)  ♦o.pqq) 

PFAO  loe*  FBACT)  *FPACT?*r)LLlM 
on  pi)  TMal ♦KATTv® 

201  »EAO  in*5,  (KAYl  (TF«IW)  #1E«1*NEXEPP) 
on  ?n4  |T«l,4on 
00  ?n4  IM«1 .NATTyP 
204  NINTFB  (  TT*  TY»)  *0 
no  206  TL«l#NLnEF 

PEAP  1P««  lOT* (TNTEP(tW) *!w«l*MATTYP) 
no  206  n««l,NATTYP 
2P4  NTNTPB  (  TOT  •  TM  )  alMTEP  <  I*** 

PEAfi  1'’'^*  .TabG(IC)  *IC«)  •NCORID) 

TPPTM  ( n  ■  ) 

00  20r  TC«2,NCnPlD 

20  0  TPOTM  (  TO  «TPnTNT(  IC-U  ♦^TAPG(1C»)  ) 

ITAPGaTPOlN'T  (^CORIO)  ♦^'TAPG  (NCOPIO)  -1 
DO  ?0?  TT«1*ITAPG 

PEAO  l-^F*  Ip  (  TT)  ,NEXtV  (  IT)  ,LAT1  tLAT?tLAT3,L0NO)  *LONG2.l0N(^3 
LAT(TT)a(LATl*  ( L  A  T  ?  ♦L  AT  3/6f) ,  )  /60,  )  BO,  017453 
LOWG(TT) « (LONG) ♦ (L0N62*L0M63/60,) /60,)«0. 01 7453 
inTalO ( TT) 

no  202  tm«i*nattyp 

SAM(TT,TW)«MNTEP(IOT*Im) 

2()?  KAY  {  TT,  Tw)  *FAY1  (KEXt  Iw) 

00  205  TCa) ,NCOPTn 

PEAO  1  OP*  CL  AT) #CLAT?,CLAT3#CL0NGi .CL0N62*Cl0N63*CLAT4 ,C|  AT5*CLAT6 
1 »CL0NG4*CL0KG5tCL0NG6 

CLAT ( TC* I ) ■ (CLATl * (CLAT2*CLAT3/60,) /60,)PO,017453 
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CL  AT  (!C.2)«(CLAT4*  {CLAT5*Cl.ATf>/#kO,»  /60  .  )  •»»  ,  0 1 7*53 
CLOKG  (  TC  •  n  «  <r.Lnss~  1  .  (CLON^J^CLONOJ/^O,  )  /An.  1  *0.0 17453 
CLC^rt  (  IC»2)  ■  (  CLO^  r,**  {CL0S65*Ct  ONGA/*n  ,  )  /  An  ,  )  *0,0  1*453 
205  CO'^TTM.5 

RFAO  ins»  (C < N) « Iw«l »^ATTYP) 
read  (PEL  (  Ii-)  t  IW«l  •'yATTYP) 

VFE»".n 

DO  21S  TC«( .KCORin 
lETWSTaTPOlKT ( !C) 
lLAST*TPlR«:T*NTARn  {  TO  -  I 

V  ALl'P  {  TO  «<',  fT 

n*AC''<!(*ilN(CLAT(TC»n)«<ilN(CLAT(IC#2n  ♦CO^  (rUAl  ( IC*  U  )  •COSCCLAT  ( IC 
It 2) )  •c^'S  <Ci.rNr,  ( Trt ’) -clong  ( IC« I  > )) 

gamma  xacOS  (  (SI  \  (CL  AT  (  iCt?)  )  -COS  (D)  •SIM  (CLaT  (lCtnn/(ST'^  (0)  •C0s< 
ICLAT ( TC»  n )  )  ) 
or  TT«TF  IWST,  TI.AST 

q.ftcrs (SIN(cLAT ( iCt  n ) •SIM (LAT ( IT) ) ♦COS(Cl  AT( ic» 1) j •cos (L AT (IT) ) • 
icns  (L'^>^^^  ( IT)  -ClOMG  ( ICt  1  )  )  ) 

ALPHA*ACnS(  (SlMd.ATdT)  )-COS{P)*SlM(rLAT(TCtn  )  )  /  (SIN  (R)  ACOS  (CLAT  ( 
I IC*  I ) )  )  ' 

RET  A«ALPHA-SA‘<MA 

TnlST ( IT) ■3**0,^ASIN(SIN(H)*SIN(BETA) ) 

DO  2n  T>««1  .MfTTYP 

IF  (S AM ( IT* 1 W) ,E0. n, )  60  TO  209 

NTMTEU(TT,Iw)3(SAP(ITtI*)»C(Iw) •PEL ( Iw) *0.5) 

GO  TC  Pin 

209  mtmtfo  (  T  T  t  Tw ) 

210  C0^TT^'lT 

212  VALUF ( TC ) ■VALUE < IC) ♦V ( IT) 

2\S  VEE*VEF*VA|.iiE{IC) 
no  221  TC«l.MCoPin 
IPIPSTsTPOTM  (  IC) 

ILAST*TF  TOST*Maoo  (  TC)  -  I 

PQXM  laO,  IC 

no  2?)  TTalFiPET, ILAST 

221  PPIM  IPS*  TD(  IT)  »V  (  IT)  tTniSTdT)  *  (KAY  (  ITtlW)  ♦IW«l*5)  *  (NIMTeP(IT* 

1  1*0  *1**1  *5) 

PRd'T  no 

nc  30n  TC*ltMCnRIO 

IF  ( ir  .fQ.NCC*^  in)  GO  TO  22^ 

PEAO  ns,  (  aP  (  V)  ,Mal  ,MAP) 

RFAH  loS*  (PIST ( I ) t I*l tNniST) 

pn  2P''  TY(*1  .mattyp 

DP  2?  ■'  m*1,kAP 

220  PEAO  loC,  {PP(I,M,IW) .I.l.vniST) 

2?3  PC  PP^  T"**!  , mattyp 
ATTCS V  (  TC* TW , I  )an,o 

222  nAWCPv(TC*IW,l)=n,p 

IF  IFJSTalPOlNT  (  IC) 
lLA‘^TaTF?R5T^MTARn(IC)-l 
IF  ( rc.EQ.MCOPTO)  GO  TO  2SS 
OP  ?7''  tt«ifi«st, tlast 
on  2?®  T*?*M01ET 
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ir  (TnT^T(lT)  .OF.DlSTd)  )  (|0  TO  125 

iuiST ilT}sl 
00  TC  230 
225  CO^TT^'lE 

lOlST (!T>«KplST 
230  CONTI^tP 

00  25''  T*i«l  ♦KATTVP 
MINCPP*NlNrPP  < IW» 

00  25''  m«2,mINC«»P 
A (M, tw) ■ (M-i ) •TNCP ( TW> 

00  2*5  TA*?.^AP 

IF  (*(v,I«) .6E,AP(I*) 1  00  TO  2*5 
NP<M,T«,)«I* 

60  TC  250 
2*5  CO^TT^'CF 

NP  (M«  i\M  )  SNAP 
250  CO^'Tt^^,'F 
255  DO  30''  T5«1  .NATTVF 
XL»'*0.0 

MT^cRP«MI'^CPP(Itt) 

DO  M«p,MlNCPP 

jRsNP (w«  TW) 

JA«JP-1 

CL'''A**'',0 

00  Pf''  IT«IfIR5T«ILAST 

IF  (  IC.NE.N/COPID)  go  to  257 

K(TT«Tto'*KAV(IT,Iw) 

60  TO  ?59 
257  KPmTOrSTdT) 

KAsivR-  t 

ppq0Pl«PlNT(0I5T(KA) « 0 1  ST (KB ) tPP <K A ♦ JA # I W) *PP C FB f JA *  I W ) * TDIST  d T ) ) 
PpOOPPaRlNT (OIST (KA)  fOIST(KB) #PP (FA* JRt I«) tPP ( FR * J0 *  I W W TO I ST  d T ) ) 
PPPOO.OTNT (AP ( JA) ,AP( JB) t^PBOBl •^PPOPp* A (M* TW) ) 

IP  (PPPCP.QT.  l.O)  PPPOB«l*« 

IP  (PPQOP,LT,').0)  PPR08*')*'' 

K (IT, Tw) «BBFAY (KAY  dT*lK) *PPROP) 

IF  (MnTEO(IT»T*<)  .EQ*0)  GO  TO  259 

IP  (PPPOB.ec.O.n)  GO  TO  259 

NU  TCP  (  TT*  Iw)  »  (MMTER  {  I  T 1 1 W )  /PPROR*n  ,5) 

259  CPTTL“(TT.Iw)«v (T T ) *K (I T t I k ) ♦•2/2 , 

250  CLi^ASsftwAXl  (CL^AXfCRITLMdTtlW)  ) 

LAMOAl *CLMAx 
LAMPAPaT.O 

27''  LAPOfi*  (L AmDA  1  ♦LAM0A2) /2, 

IF  (  (1.  AvDAl -LAMDA2)  .LT.O.nnl )  60  TO  ?R3 

call  FnAMdh) 

IF  {  ATTACK-A  (t' t  IK  )  )  275*267*260 
275  LAMDA1«lAM0A 
60  TO  27') 

260  LAMDA2«LAM0A 
60  TO  270 

283  LAmOAwLAMDAI 
CALL  FCAMd^) 

ATFMPaATTACK 
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DTEMP«nAMA6p 

AMOA? 

CALL  rnAM(Tto) 

AlaAPS ( ATEMP-A (M, iw) ) 

A?«ABS(ATTACK-A{m,IW) ) 

If  <*?.Lf  .Al  )  (^n  TO  385 

ATTAC«»ATEMp 

OAMAQPanTEMp 

285  MCaM.l 

APTaA (8t  In) •IKCB ( Iw) 

IP  < attack, fiT,  aTTCRV(TC»IW|MO» •ANO.ATTACK, I  T, API)  00  TO  ?87 
ATTCBV (lC,Hi,M,,ATTCRV{IC.IW#MO) 

OAMCBV  (  TC  •  Tm*<)  aPAMCPV  (  IC  •  1P«M0) 

00  TC  ?BP 

287  ATTCRV  (  tC«  TP»**)  aATTACK 
OAMCBV ( TC# TP fM) an AMA6F 

288  CONTTnLF 

IP  ( ATTACK, FO.O. )  GO  TO  2Qo 
XL**?«nAMAGP/ATTACK 

XL^aAMA* 1  ( XLH#  XLM?) 

290  CO^TTMLP 

291  LaM^^A*  {  TC#  th)  axLM 

IP  UPBTNT (1 ) ,FQ,0)  00  TO  300 
MINCRPaMNCPP  (  IW) 

00  298  Mal,pINCPP 
X(M)aATTCRV(IC,!P#M) 

298  V  (M) aOAMCBV ( IC# lP#M) 
print  n5#  TC#TP 

VMTNan , 0 
VMAXaV (MisCPP) 

N3a  \ 

LIKPa-» 

CALL  PLOTn  (Y#X#YMIN,YMAX,N3#MINCPP#LT^f ) 

300  CONTINLP 

CALL  ^prONO  (TIMP) 

PRTNT  190,  TIMP 

301  RpAO  105#  ( WO ( TW) # IKa 1 #NATTYP) 

IP  (pop, 5)  999,302 

30?  00  303  TWal,NATTYP 

IP  (wn f Ta) .LE.NWPPMX ( IW) )  GO  TO  303 
PRINT  199,  IW#NWPPMX ( IW) ,W0 { IW) 

GO  TC  3ol 

303  CONTINUE 

no  304  Taa| ,natTYP 
L^iMHCA  ( iw)  *n,0 
00  3'’%  ICal,NCORIl) 

304  LAMHPa ( IW) aAMAXl (LAMBDA { IW) # L AMMA X ( T C # I W ) ) 

I P 1.  A  G  »  ^ 

CALL  AlhCAT  «NATTYP,NCORIn,NINCRP# IPLAG) 

0X8»n , no 

no  315  TW8l,NATTYP 

IP  ( TOTAL* ( IW) ,LT,wn ( IW) )  QQ  TO  305 
DL AMOA ( Tw) aO,05 


m 


•'t 

4W 


■si 


41 


30S 

3in 

311 


3an 


3?? 


3?3 

3?5 


330 

33? 


3j5 


GC  TC  tin 
OLAMOA (!•) ■-o.nl 

AMRn4(lM)*(l,*nLAM0AltW) 

CONT  T^^  r 
ITru.i 
Tl* J«|  <',F4n 

CALL  AL^CAT  (NaTTVP.NCOPID.I^INCBP,  IPPlNTf?)  » 

CALL  orsTcv  (\aTtvp,ncOPIO) 

ir  { TOPTNT O) ,pQ,n)  flO  TO  3?3 

AT^TAi bO.o 

OTOTAl an.o 

LT0T41 .n.o 

PPT^-T  MS#  TT€p 

PPINT  i?n 

no  3??  T^»«1,^ATTVP 
AT0TA(  ■ATDTAL^TOTALA ( !W ) 

0TnTAL«nTOTAL*T0TALn(I»l) 

LTOTAt«LTOTAL*TOTALLni) 

I  TOTAL!  I  TW)  •lambda  (  tW)  tWD(  fW)  .TOTALA  (JM)  »TOTaLD  f  JW)  , 

PPINT  130.  aTOTAi  .OTOTALtLTOTALtVFjr 

nLMAK«»*,0 

DO  3?«  TM»1,NATTvP 

0LMAX«am4¥1  (OL»<Ax.aBS(DLAMOA(Im)  )  ) 

IF  (OLMAX.LT.DLL IM)  GO  TO  A50 
3tn  tw«i,natTvp 

(TOTALA ( Tto » .NF.wn ( IM) )  GO  TO  5no 
TC  ?n5 

3AR  TM«t,NATTYP 
(DLAmDA ( IM) ,GT. n, )  60  TO  3A5 
(T'^TALA  (IM)  ,GF.mO{  IM)  )  GO  TO  340 


DO 

If 

on 

IF 

60 

00 

IF 


3A0 


345 


350 

355 

360 

365 


IF  .  _ 

DL  AMO  A  ( TM)  «PPA(:T)  *r)LAM6A  (  IW) 

60  TC  34O 
nLA^'OA(TM)«FPArT?*nLAvr)A(Iw) 

60  TC  355 

IF  tTOTALA ( IM) ,GT.MD( IM) )  GO  TO  350 
OL AMP A ( IM) ■FPACT?«OLAmDA I IW) 

GO  TC  3*0 

Ola MCA ( TM)«pPArTi*nLAMOA (Im) 

Gn*??^:SI?’ **''**’  frii-^A*/2o..AMiNi{nLAMnA{iw),o.5)) 

‘-'^L^A><''20..AMAX1(DLAM0A{tm)*-0.B)  ) 


TS»u  •  n 

TA.^.n 

no  37'^  TM«1,NATTyP 
Ts»Tt*i  AMPOA ( IM) *T0T4lA(  iM) 
370  TA»TA#L AMPDA (  IM) •MO  ;  IM) 

IF  (CX,6T,o,)  oO  TO  3P0 
IF  (TP.TT.TA)  no  TO  375 
0X*FBACT 1 *0* 


00  Tr 

3?^  0x»FOAr 

A  A  T  A  'Ion 

300  if  IVs.ot.ta)  r,o  TO  305 
n<«F0AC  T2*n« 

00  TC  ■»Q5 
300  ni»FBAfTi*n» 

3«P  n*«AVA)l  1  tOLMA»/?n,  •AMlNUnitO.S)) 
on  Tr  400 

395  D*«amtm  (-r)L^**''?o  ,«A“A»l(nx#-0*S)) 

♦  or  on  41''  tw«i,natTvp 

410  LAMPOA  (  Tw)  •(  1  .*01  AMOA  )  •  { I  .  *0*  1  *1  A“BOA  f  Hi» ) 
ITF0«TTF«*1 

ir  ( ITrq.GT, I TFOPX >  GO  TO  450 
GO  TO  3?0 

*«in  no  AOf'  t4»i*^attyp 

LAMprA(I>*)«LANnAB(T>i) 

400  TOTALA  (  Ti^l  ■TOTaG  f  T  <() 

00  4?^'  TC«1  ♦NCOBjn 
ITVPe  f lO ■ITTPF«< IC) 

PtoEAP ( TO) •PWEAPG  r IC) 

FOAP  (TO  ■FOAPG  (TO 
470  TNTEBV  (  TC)  ■TMbvp  (  IC) 

ITFP.TTFQB 

500  call  close  (NATTyP,NCOPI0tNlNCOP) 

505  00  T“«N'^‘‘TTvF 

TOTALO  (  TIk  )  «o  .0 
51  0  TOTAL!,  (  T«i )  ■0#(’ 

no  5?f  TC»l •NCORin 

I  *  1 7  V  p  P  ( I  r  j 

IF  (T^.PQ.'')  60  TO  5?0 

TOTALr'(Ti«)«TOTALn(I)i()*EnAM(IC) 

TOTALL  (IW)«70T4I.Lnw)*E0AP(  IC)  -LAMROA  (  IW)  aBWEAP  (  IC) 

52P  continue 

PBTnT  140#tTEP 
ATOTAL«0.0 
nTOTAL*''.  0 
no  541  T«<*)  ,N4TTvp 

total  A ( IX 1 aO . 0 

54l  TOTALO  (  IX)«0.'J 

no  ftO''  tc«i*ncortp 

Tv<«  ITYPP  (  lo 
IF  dw.FO.'')  60  TO  59R 

XEAPn*R'*EAP  (  IC) 

CLMAXbO, o 

IFIRST«IPOIkT ( IC) 

ILA5T«IFIRFT*NT4RG(  lO-l 
IF  (  ir.FO.NCORin)  60  TO  547 
DO  S4?  IA«P,NAm 
IP  (XEAon,GE.AP( lA) )  GO  TO  54? 

NPn*TA 
GO  TO  543 


54?  COM!M,» 

NPr>aS  AC 

543  CC\TlALr 
JnaNPri 

547  00  546  TT.TC  1««;T,1L4ST 

If  (  !«*  ,PQ.A.C0«in)  00  TO  544 
KR«iri4t (  lY) 

KAaKC.I 

pppom .OIM (0157 (K4) tOlST iKB) .PPtKA, JAtlW) tPPfKfl# JA*1W) #Tnl5T(lT) ) 
PPROPP.BIAjT  (0157  («  A)  *0 1 57  ( KB )  tPP  (K  A  *  JP*  I W)  tPP  (K«  t  JB#  I W)  •  TOI ST  ( 1 7  )  ) 
PPR0P«CTN7(AP( JA) .APC JBI (PPPOBl tPPROB?»MFAPO) 

IP  (PPBOR.07.1 .9)  PPROflal.O 
IF  (PPBOR,l7.0,0»  PPROHaO.ft 
K  (  IT,  T*>  aBt-KAV  (KAY  (  TT,IK)  ,PPR0B) 

IF  {MKTEB(T7»TK)  ,FO.O)  00  TO  544 
ir  (PPB0H,F(5,0,O5  OO  TO  544 
MAjTFB(TTtTw»«(AiP  TEPUT*  lK» /PPP0B*0,5) 

544  COITi  »'(T7,Ik)«v(IT)«K(IT,Ik>*c2/2, 

CLMAXaAMAYl (CL“AY,CPITLM(IT,IK) ) 

546  C^'^TINlE 

LAWOA 1 aOLKA Y 
LAMnA?a9.0 

545  LAMOAa (I  AMOAl *1 AmHA?) /?, 

CALL  POAM  (fW) 

IF  (  aTTACk-kEAPO)  550*5BR*555 
55^^  LA^OA^aLAMDA 
(50  TO  46^ 

555  LAMOApal.  AMpA 

569  IF  (  (la^OAI-LAmOAP)  .LT.O.Oni )  GO  TO  579 
GO  TO  445 

579  CALL  CLOSE?  (latWEAPO) 

580  PBIN7  160,  TC»TK 

PRINT  162,  attack , damage, value ( IC) 

ATCTAL«AT07AL*ATTACK 
OTOTai apTOTAL^nAMAGF 
TnTALA(TK)«T0TALA(TW»4ATTACK 

totald { IK) ktctalp ( iKi ♦damage 

EDAM ( ir ) kOamaGE 
00  50''  TTbIF  IPST,  ILAST 
550  PRINT  165,  TO (  IT)  , ATT ( IT) ,nAM( IT) 

GO  TP  6nG 
595  PRINT  170.  TC 
609  CONTTNLF 
PPINT  1<54 
PRINT  ISO 

DO  603  ICal ,NCPPir 

603  PRINT  155,  TCf TTypf ( IC) .RkfaP ( IC) tEnAMl  IC) .value { TO 
PRINT  193 

no  6N4  IMSI ,NATTYP 

605  PRINT  195.  Tk.TOTALA ( IK) .TOTALP ( IM) 

PRINT  196.  ATOTAl .DTOTAl , VFF 
GO  TC  39l 


pbtnt  jq* 
ion  ro3«*T 
in§  rnpMAT 

10^  Fco»**T  (?iin#Fi'>, 0,2(3*, II, ?i?)  .Srs.nj 

10*  F(^OmAT  (♦(  3*  #F  9F2,n  )  ) 

1 1  r»  FoowaT  { 1-1 ) 

115  F^OMar  (iNt,*  itfoaTTOn  no»<BF»  ••!?> 

l?fl  FCOwiT  (//•  n»  r»LA*<«nA  WEaP  ofs  WFAP  OBT 

1  paMACiF  uanwAN*/) 

l?^  FCPMftT  (In  ,(tft,A,2Fin,%,?rlP,ft«?Fl?,4) 

130  FCP'^AT  (///•  Total  ATTACkFO^  ■  •#F4,n/*  total  OAPAOE  ■  •*F10#4/ 

1*  TOTAI  LAOPASO,]  m  •»F10,A/«  TOTAL  VAH'F  ■  •*^10,2) 

134  FfjOMAT  (////•  ropoToOP  NijMpfp  ••!?/•  wEAPom  MiwpEP  BtlP) 

|40  FCPmAT  (lHl,*FTrjAL  ALL^'CATTON  FOUnjO  AFTE«  n.0SIN6  ON  ITCPATION  •* 
112) 

150  F0P«AT  (/////•  CnwoiDOR  HEAP  Typ  ATTACK  DAMABf  VALUE* 

1/) 

154  FOPi-AT  (In  ,  I  5 ,  1 1  2  ,  F  I  ?.  0  ,  F  1  ?,  ?,F  I «  ,  0 ) 

l#,0  fopmat  (1ni,«fikai.  allocation  to  TapgfTS  TN  COBPIOOP  •*!?/•  This  c 

lOPPlcrn  «AS  ATTACKEO  RY  wEaPON  •tl?//) 

format  (•  total  ATTAChEOS  to  rOPPinoP  ■  •*F4,n/*  total  damage  to  C 
lOBPloOQ  a  a.Fl^.A/a  TOTAL  VALUE  OF  COPPTOOP  ■  ••Fl0,2//*  TARGET 

2  ATTACK  OAMAOE*//) 

165  format  (In  , I S , F 1 2 , n , p 1 R , 4 i 

170  format  (  lNl  ,4C'^BRI00P  HAS  NOT  ATTACKED*) 

mo  format  ( Ini ,«COWQinOB  NUmreR  *#12//*  in  value  dist 

1  K**4n*«*TK TEbCFPTORS*//) 

185  format  (In  ,  I5,?f1  f),2f5Flo,4,51«) 
iqfl  format  (•  time  a  aiFlr.A) 

WS  format  (////•  WEAR  Typ  attack  OAMaGF*/) 

iq4  format  (IMl,*  summary  OF  FINAL  ALLOCATION*) 
lq4  format  (1m  ,I5,F1?,0,f13,.'  ) 

196  format  (///•  TrtTAL  ATTACKERS  •  **E5,n/*  TOTAL  DAMAGE  •  *#E12*A/ 

!•  total  VALlE  ■  •#F12,4) 

19R  format  (1-1, •error  in  Input  values,*/*  number  of  aTTackEP  types  e* 

ICEKO*  TmF  nlmRfp  of  CORRIDORS*) 

iqq  Format  (1H1,*E0R''0  1^  INPlJT,#//*  TMIS  CASF  WILL  ^OT  BE  RUN,*// 

!•  mjvhfp  of  weapons  of  TVPF  ♦,Ii»#  cannot  exceed  *#I5// 

2*  This  case  REOUFSTFO  •tF4,«)*«  WEAPONS  OE  TMAT  TYPE*) 

999  STOP 
ENO 
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c 
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F 

A  I  I  '  ■‘I  •  • 

I  M/.  4.<  S')  ,  '• 

i  1  J  ^  >  ;  I  »  *  I  ♦-  I  W  S  T  ,  I  I,  4  ^  r 

f  ^  (  L  *  M>)/  •  I ’►  •  C  <  I  T L**  <  T  f  ♦  I  •  '  '  y^^,l 

A  f  r  «  u  » 1 1 . 1  )  •fti  '.  ( y  ( I  r i  “K  ( IT.  j 4) ••?/ ) )  ••^ 

4  >  T I  « I N  f  ( •»  f  r  1 

1 1  .  M  I  a  I  <Ltv(lM.'^(!l»T»i.'TfJ) 

1  r  (  ■<  I  4  I  fc  •<  •'  I  T  ,  1  *■  »  ,  w  >  ?<?  » 

I  n  «i )  A  4  I  /  (  4  f  I  I  •  ^.  I ' .  f  t  ^  (  I  f  ♦  I  •* )  > 

I»-  (I  I  ,L.‘ >’  iA)  jv  f''l 

??  >  4  f  T ^».  r  r  1  •  I  . 

I  ,  a  4C  t  »  (  f  ’  )  T  f  •  I'J  ♦  ‘TT/*) 

!  ^  !  a ,  1 A  4  j  « L  ••  <  '  *  A  r  T  I 

hii'-AI  f'' 

ii  fC  ^4' 

A  f  T  I  I  T)  aft  f  M  •  4l  a  (  T  T  «  !•' ) 

i>*\ (  f  T  )  af'  k  M 

...  1  1 

P4.’  4 1 T I  f  ? )  «/■  r  r-^ * , !  Tr  .M  T  r ,  t  o 

O  (  I  I  )  af  ^ 

PA"  A  r  r-kCKaA  1  r.-CA  •  A  f  T  (  T  f  1 

i'».  at)  A  *..  iw  *  >  '  <  (  M 

' .  J  I  I )  )  0  f' 

?‘)0  A  I  T  M  r  )  at) .  I 

I  It  *  k  {  f  r  )  ai> « 

I  1  "kn  C  •  I 

s  i>.  r  t 

(f  I  - 

I' 

j: 

is;' 

[  ► .  '\'*  r  I  vj\  .1  I'-''  ..C  {  .  ♦  A 1 C  ■ 

f  i  (Ta  A*  (  ,  ;  I  .  .tas  i-' MC  )  )  •r  <  •’ (C  )  )  ) 

I  w .  T I  •  ' 

i.  ^  . 

I' 


I 

% 

I 

# 


i 


t 

')■ 

A 

■i 


ii 


I 

i’tt' 

■« 


f  )K.t n-M  '♦HtSifC  'KuV**-’-^) 

I  )«pi^*  ( I  •  ■'*'  A*^!*  Ka  f )  «  ( I  •  •>  •  ii^Ar)) 

»»U  ' 

I  fs  I 

•%  T>  4k>a(^  4P  I  >  «  ) 

I )  I  «  »  .  ■’-  Tt "  t  '  •  ♦  '  )  “'^ 

!►  (iiSt  .'l>  .If,  l,u^-5)  r,i'  T.) 

I  ►  'If  ,  'j !  .  ■>  ’ )  t’  I  T'j  (?»• 

! I  * :  r* » 

\  I  /nt> 

'  •.  I  '  )  -* 

1  I  r  <  1  y  a  a 
»w  ►_  I  I  »<  ■ i 

?.i  M'*?  .r  ii)i>  «a«  >«i  l 

T  t-^  i 

►  •  ' 
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►  JN"  r  I  ‘)rj  ►  'il  (  ai  »<?»Yi  ♦V?«a) 

^  prok  jrl^S  L''iK  I  >% 

-Jl  M  »Y  I  -  ( -V  I- a)  •  ( Y^-y  I )  /  (  x?-y  I ) 

-  T  i  I W  ', 

6-  . 


! 


1 


c 

c 


1 


■» ' 


1 

1^4 


«;  I  ■•>'  '  j  I  IN*  P..ir  (  y  1  ♦  An  «  Y*^  I '•  YMAXtl  J»NP1’ *1  IN§.  ) 

y  1  ♦  \  y  na  ’tS 

NP'J  a  NO'-’  >i  -y  yjf’  POJnT^  *^L'''TYf-.n 

r, !  >><►;  )  p.  J  I T  i lo  1.) 

I'll  <►.  ^sr>'i  1  ( i )  •  A*'  ( I ) « -'T  Ic^  n  1 ) 

I  f!  »r  r<  S>'  AC^M  '  P  ♦  A'»Y 

i'ATA  <l_aN^  »  '  lO  f  »  /  I '1  «l-',«lw*/ 

!►  (  NMP.t'J.  A  )  '^*  r«)-<N 

Y'aA),/(Y'lMN-r^;N) 

Y  (  KK  (  I  )  rYMT  4 

Y  >1  .A»i  n,/y. 
lA' '  A  'A  I  a  I  »  S 

Y  I  I  r^  (  I  ♦  U  =v  T  ;  r;Y  (  I  )  ♦  Y-A !  ., 

I  I  r  10  ** .  {  y  U  (  n  •  I"*  1  ♦M 

►  P^AT  ACK*,F?,  0  10,0) 

wi^Tfr  '  A-*) 

pr  M  -A  ?  (^-»<'-lIY(l  . . ■)  ) 

NLS  '•  J*  yPP 
L  I  *l  I  ♦  I 

IA.I4  I  a  1  ,'5  J 
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J 


c 


13 

13-* 

1? 

n 

i') 


i 

\n^ 

4^ 

5'' 


(  I  )  m'il  .\ 

S*-  »  r  I  'i  *  !  i  L  1  <tS 

J  »  ’  (  I  ’  f  I  I  •  I '? «  M )  t  L  ; 

'  '  M  )  I  «  1  I  »  1 1  «  I 
S^4':it  {  I  )  »  )  )T 
( ■>»)  *i  oT 
S'^ATC  ( I ) ».)  ' r 

>4iJa  \  ,  Ml.  «  n 

H«/>w  (  jJ) 

Mr  1  a  I 

I  M  =  «  Y  I  (  -T  M  I)  •Yl  *  U  ■> 

!»•  n  T  I  .LI-  ,  ’  Y  I  .bc  ,3^) 

!►  (*''  I  J,  )  )  T  '  1 

f  »^w^Cr  (  I  Y>  ) 

(  lYII 

T  fP  (  O*  i  H)  r  I  (  MkJ »  ,  X  «  !  S*^  Cf  (  T  )  • 
Fi-VM4T(^<rl3.*.fn,.),UsUn 
IFImCI.mK,  -)  'Cf:  (  1  Y  I  )  alSJl 

C  J  iri  ^  it 

T  r  t-  ( f> ,  I  '  ^  1 
►y  •  T  1 1 W  ■ , 

F  > 


1  ,Sl  ) 


c 

c 

c 

c 


S  H-»''ijT  1  alJC  »r  (\4T  T  YP.NjrOW  IO*MNCnP»  IFLAti) 

CUM  -«o\/CAi.L  UC/  (^)  « I  type  (15)  .WiKEAPi  15)  tEOA^  ( 15)  * 

Ml  TrPw(l'i»S»^l)  .DAmCPv'  tLAMMAX  (M«5)*lMTEPV(l3)» 

^T"T  ALm  (‘^J  «  I  )  r  \i.  MS)  *  TO  TALL  (5)  «  VALUE  (15) 

UME  'ISIU  *  f't  ^C.P'MS) 


Hr  Al  l  Am  <  ) » «  L  a  '^■iA  X 

iA  1  (s I *''iA rr »p 

TuTmLA  (  I)  =". ») 
ruTALJ ( I ) *C.O 
T  jT  .LL  (  F  )  siA.O 
nu  M)l  IC*  I  fMC  MIO 
T  r  YUf-  (  iC  )  a  ' 

P*F  ■.H(  IC)  *1’),') 

EuA  M IC) *  J,  0 

TmTFWv  (  IC) 

Y-  AYSaO.  A 

fm  p-^r  [M3 1  ,,NA  r  fYP 

[)•  (  LA -(SHd  ;  I  4  )  ,  AT  .L  AMMax  (  I  c  «  I  W)  )  GO  TO  ?hi) 

ULA  M  30,0 

M  I  «jrHPxN  {  f  I) 

I  >'.>  710  K37  ♦  H  1 

HU  AiixUA'^CHV  (  IC  •  I«»  ,x)  -L4M-)r  A  '  I  M)  •AT  TCWV  f  TC  »  1  «»|^  ) 
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!►  t-^LAi.L- .-^L' »  G-’  >!■' 

<>  '•  I  *HL  ^  I 
><«-•/ 

?  I  '  r  >n1  I  MUt 

f  rC*<v  (  tCt  I  •»<^^) 

W'  i»i'  ^'1(:w V  ( I  r  1 1 «( ♦  AK ) 

If' 

’4  .  S  >1  I  (  t  4  I 
VI  ?yo  i*sP 

!*•  (  ,  'TTc-vv  (  IC,  )  uO  rO  ??0 

.1  It 

??  t  <:  I  |T  1  ,nf 

t  1^ .  1  ^  -  I 

WJ«<INI  fAI  !C»1»I»^J)  ♦ATrCf*V(IC»l'<«KK)  lOAHC^V  (  ICf  lot*^  J)  • 

t  1  \  >4C-<W  (  •  I  4  •  I  •**  \  ) 

•^L  ^  ^  I  ■^' i\il4|*WA 

?s>'  !►  (  ^L.A4i  ,i.r  •  wuXS)  GO  T  1  r  » 

T  I  t'*»  (I  ' »  *  I  << 

P  4P  I  TC  I  ■»<  A 

‘  ^  M  IC  )  a*^  I  — 

f  iT  4^  (  I  “  )  aA^ 

^  •  A  «Sa4L  4  ^  I 

> 3  '  f  t  ^ilP 

?  •  a  ^  f  »»*f  »  I  ) 

f  Ta.  i*(I*»*TGT4i.AIIaJ*>4  4FflU(lCl 
▼  f\C  rA:0(I«»*fc  )AvnC) 

T  »A^^(!*|aT  '  (  At  L  n«  '  ♦  S 

1  '  *  '  *1  »  I  s.jg 

I>  i(fuA««;w»  '  Gl  Tl 

»  .1  i  J*! 

M  '  ic* '  •  •:  •  <i  ' 

^1'  h>«isf  il^*  tw»i‘'r*^F(I'!«’**^A*^ClC)*EOAM(IC) 

^<'l►  •}<A^  fT*'*’  4Eap  oamage*/) 

t  41AI  '«’:>» 

T  1-3  V 

c  . 


<;  I  .r  "t  >T;,t  (nat  rr»^«  c  'Hinj 

c 

r  »  •*-«')  t/C  M  l- '<- /<•  >  t-i)  M\  (S)  •  I  TYfc»£  (  Ibl  .rtgitA**  <  •R  IS)  » 

t  A  t  rr.j  ^  (  »  S»  ■'»<  I  »  *  )A»4C<»  (  I  ■>,'>1?  1  I  (lo^SJilNTE-^VClSI* 

■>T  »T  At  A  ( s)  «  »  it  AL  M-I)  »T'!!  ALl  JS)  ,VAI,  »E(1S) 

L  1 4  HI  j/c-t  s  r  / » «  )•  T  jf  A-u  3) . !  rt^»Es« )  o)  iMdEAPft  IS)  tti'A'^^ns) . 

I  T  v.T-<  VM  (  \  s;  «L  A  »  i  ^  (  S)  *  I  1^  f 

C 

Pr-  ftl.  LA**'*'*  '  *1- A**  1A-) 

c 

c 

T-'i  J1  *'i,0 

lit  >  J  J  *  I  *  ^  A  t  ^  Y  ^ 

!►  :  . j(  j  j)  .f  j.  '.0)  I'O  fo  '»u‘* 

f*  ii«l4jl*  (  (  r.iT  iL«(JJ)  -  f'Ai)  t)  )/«r)(jj)  )  ••i 

t  r  )  <>s  ^ 

?-'  )  T*  Jl  a  r  3  Jl  •  T  ' T  \L  \  (  JJ) 

-t  T'  I'l  f  I 

I  )•■  «  ?  a  J1  .  >  r  ,  T  »  J  t  ''.0  T  )  -i  , 

I  Ik  !a  I  •  J I 

'  '  ■»>»  )  j  Is  I »  .  \  1  r  yo 

I.  .i-^'AfU  JJ)  sLAt-*  )A  (  JJ) 

<k'V)  r  'T  •*  S  (  J  J)  S  I  I T  \L  ■»  (  .1  J) 

In'  »*  T'  i  J  a  I  »  >*r' '  H  1 1' 

I  I  yPE  u  JJ)  *l  f  :  (  )  J) 

AP  ■•  (  ,J  J  »  sp-  -  A  J  (  J  J) 

F'.'A  (  JJ)  sKil  A  -i  (  IJ) 

*«;a  I  -<  (  j  J)  si  ifF  W  ijj) 

I  (  ^  « a  1  T  > 

S  '  '  f  -  Ti  M  -) 

F  ' 
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c 

Cu-4*  M/C*Lt  >€/•  >(>)  vm  (S)  *1  TrP*: « |S)  •Hw^aPC  is)  tt  • 

MrTC^V(lS#'5«<>l)  tOAMC^Vl  IS,S.?|)  tL4MMA>l(  130)  tlN’EUVIlS)  f 
?TuT  \LAO)  .  r  )Tal  MS)  «T')TAH  (‘<')  *)» AL'»6  <  1  S) 

c 

tM  ■<t  '•^lOS 

c 

w»>  flL  LAM-»i)4 
C 

c 

0.1  nsNNATTyP 

( roTALA(  iw)  ,uH:.wo(i<()  )  1,0  to  300 

iV.CaTUTAt,**  (  1  <«  )  -  <'M  TW  ) 

SUO^tslO.hAr, 

no  ^SO  IC*i  ♦■xro  ^TD 

IF  (  ITYP«i(  IC)  .  M  'iO  ro  >*‘f' 

?io  IF  ( iNTE-^v  ( ic)  .FQ.  I )  (’tO  rn 
KK,IHIEPV( ID 
t^jsKK-  I 

IF  ( -VTIC  <V  { ic«  I -<»KK)  .Nc.fl  ITC^^V  ( IC*  I>^»KJ)  )  00  10  2?n 
TmIF*.(V  (  IC)  sI  MTF.MV  {  IC)  -1 
T  )  21*' 

??■»  SLOPElx  (D/',MCHV  (  IC*  I«»»KD  -04MCRV  (IC*  IW,KJ)  )  /  (  aITC«V  (  lC*Irt*KK)- 
1  ^  I  Tr».W  (  IC*  I*  *^  J)  ) 

IF  (SL!»Pr.  I  .r.  r,SLOF>F)  bO  To 

TC.)tC-<*IC 

CvjnrT\'jE 

ICsTCUtC-^ 

^^*1  'iltpy  (  10 

K  J  SKK  «  1 

OlF-r  aUxEAP  (  IC)  -  M  I  \/  (I  C  *  I  *  K.J ) 

IF  (OIFF.-iE.UECl  bO  10 
7  »T  41.  A  (  I  *)  sTO  (  <\t.A  (  ivn )  -,jtF  F 
••'»f:ap{  IC)  x.inc-'/  (IC*  I'^jKj) 

F*)A-M  IC)  aiUMCOv  (  IC*  I  V*SJ) 

IF'  (0^tAP(  1C)  .F  J.r.n)  ITY-^F  ( Ic) ®o 
I*.  T-*  (  IC)  *  INI*-  ov  (  TO  -1 

f'hrsOtC-t)!*-^ 

AO  ro  20 '7 

?<,o  w*>FAP(IC)='<''EaP(IC)-()EC 

EC  4  K  IC)  a»^If  T  (  aTTCRV  (  IC*  I  'M^^J)  ♦  ATTCHV  (  TC*  IW*KK)  ,DAMC^V  (  IC*  IW.Kj)  * 
|i)A  ^ChY(IC*I*’*<K)  ,-»^EAJ(IC)  ) 

IoTaLA ( IW)  =  r JT ALA (  IW) -JtC 
3m  COOT  I 'JOE 

DO  *00  I*»=l.  OATIYP 
MIiNlCrJP*Nl  mC'^P  (  I  ^) 
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